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ABSTRACT 

The excess numbers of blue galaxies at faint magnitudes are a long-standing cosmo- 
logical puzzle. We present new number-magnitude counts as a function of galactic 
morphology from the first deep fields of the Cycle 4 Hubble Space Telescope Medium 
Deep Survey project. From a sample of 301 galaxies we define counts for elliptical, spi- 
ral and irregular/peculiar galaxies to / = 22. We find two principal results. Firstly the 
elliptical and spiral galaxy counts both follow the predictions of high-normalisation 
no-evolution models at all magnitudes, indicating that regular Hubble types evolve 
only slowly to z ~ 0.5. Secondly we find that irregular/peculiar galaxies, including 
ON multiple-peaked, possibly merging, objects, have a very steep number-magnitude re- 

5-H lation and greatly exceed predictions based on proportions in local surveys. These 

^ systems make up half the total counts by / = 22 and imply the rapidly-evolving 

component of the faint galaxy population has been identified. 
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INTRODUCTION 

o 

_ Counting galaxies as a function of magnitude is a tradi- 
{■C) tional yet active tool in extragalactic astronomy. Originally 
C ~_1 a cosmological test (Hubble 1926), it is now recognised to be 
more sensitive to galaxy evolution (Brown & Tinsley 1974) 

ON 

for which evidence is accumulating. In particular, counts of 
Q galaxies in the optical passbands have revealed a dramatic 
O | excess of 'faint blue galaxies' (for a review see Ellis 1990). 

I Extensive redshift surveys have been undertaken of these 
£~ objects selected in B (Broadhurst et al. 1988, Colless et al. 
^— ' 1990, Glazebrook et al. 1995 A, Ellis et al. 1994), / (Lilly 
—j 1993, Tresse et al. 1994) and K bands (Glazebrook et al. 
1995B, Cowie et al. 1994) which have been used to construct 
the respective luminosity functions as they evolve with red- 
shift. At short wavelengths, where the evolution implied 
by the counts is strongest, there appears to be an increase 
in the space-density of galaxies at luminosity Mb ~ —20 
(we use Ho = 100 km s" 1 Mpc" 1 ), over < z < 0.5 (Ellis 
et al. 1994). At longer / and K wavelengths, the excess ap- 
pears to occupy fainter portions of the luminosity function 
(Glazebrook et al. 1995B, Cowie et al. 1994). 

There are still some major uncertainties of interpreta- 
tion. Firstly there is the fundamental problem of </>*, the 
local normalisation of the Schechter luminosity function, 
which remains uncertain to within a factor of two. Clearly 
the higher the local galaxy density, the less evolution has 
to be invoked. Local surveys have given a low value; Efs- 
tathiou et al. (1988) and Loveday et al. (1992A) both find 
<f)* ~ 0.015ft 3 Mpc -3 based on redshift surveys to B = 17 
where the median redshift is z ~ 0.05 and the effects of 
large-scale structure should be negligible. 



Adopting this normalisation leads to a factor of two 
excess of galaxies as bright as B = 19, where the median 
redshift is only 0.1, increasing to a factor of four by B = 24 
and z = 0.5. This has led Maddox et al. (1990) to ar- 
gue for extreme evolution at low-redshift. On the other 
hand Metcalfe et al. (1991, 1994A, 1994B) and Koo et al. 
(1993) have argued that the bright data is not representa- 
tive: possibly being subject to selection effects, large-scale 
structure or problems with the local galaxy luminosity func- 
tion. They prefer to normalise at B = 19, equivalent to 
</>* ~ 0.03ft 3 Mpc" 3 , which still leaves an excess at B = 24 
of x2 which requires evolutionary explanations; it is not pos- 
sible to explain the whole of it via local effects. Regardless 
of </>* the slope of the counts is steeper than the no-evolution 
model over 15 < B < 19 and this points to some anomalous 
behaviour in the local data. 

Secondly the excess at B = 24 is now established as be- 
ing caused by an increase by z = 0.5 in the density of objects 
with M B = -20 (Glazebrook et al. 1994, Ellis et al. 1994), 
i.e. of similar luminosity to local giant ellipticals and spirals. 
What has caused this population to disappear by z = 0? 
Proposed explanations include: (i) extreme differential lu- 
minosity evolution in which dwarf galaxies are preferentially 
promoted from the faint end of the local luminosity function 
(Broadhurst et al. 1988) at z ~ 0.5; (ii) merger models in 
which the extra population is merged into the local popula- 
tion by z = (Rocca-Volmerange & Guideroni 1990, Broad- 
hurst et al. 1992 (BEG)) and (iii) 'exploding dwarf models 
(Cowie et al. 1991) in which the blue galaxies are a new 
population of small starbursting systems: the starburst is 
powerful and dominated by high-mass stars; the remaining 
gas is expelled by supernova driven winds preventing further 



2 K. Glazebrook et al. 



star-formation and turning them into dark systems by z = 
(Babul & Rees 1992). 

Ideally we would like to directly identify the population 
of objects contributing to the excess and track its evolution 
with redshift. This has been attempted using the strength 
of the [Oil] line as a star-formation indicator (BEG) and the 
population of such objects in the B luminosity function has 
been found to evolve rapidly (Ellis et al. 1994). However 
[Oil] emission may be a transient feature and it is more 
sound to examine structural morphology — this requires 
sub-arcsecond resolution to see kpc scale structures at z ~ 
0.5. Limited work has been done from the ground in 0.5-0.8" 
seeing; this shows that a high fraction of the star-forming 
galaxies to B = 22.5 have multiple cores (Colless et al. 1993). 

OBSERVATIONS & DISCUSSION 

The new data we describe were taken as part of the Medium 
Deep Survey, a key project designed to maximize HST sci- 
ence by taking WFPC2 or FOC frames in parallel, while an- 
other instrument observes a primary target (Griffiths et al. 
1994 A). Thus a large sample of random fields results with 
exposures of several thousand seconds and reaching depths 
of / > 20 (note unless stated we use / = F814W). Our 
particular sample consists of all 13 deep fields with total ex- 
posures > 4000 sees observed with the WFPC2 camera since 
the start of Cycle 4, when the spherical aberration was cor- 
rected. The data were reduced in a standard way: pipeline 
bias/dark/flat corrections were applied and the data aligned 
and coadded using a cosmic ray filter. Cosmic ray removal 
is especially critical due to the high incidence in HST data; 
we use the 'crreject' algorithm as described in Griffiths et al. 
which is found to work well for WFPC2 images and our sam- 
ple excludes fields with only 1 or 2 individual exposures and 
galaxies on the PC chip which has a different pixel scale. 
The resulting total area is 51.3arcmin 2 . 

Object selection is another critical issue for HST data, 
as a galaxy which produces a single image from the ground 
can when resolved, have a complex structure. For example it 
would be incorrect to split spiral arms into multiple objects 
and it is unclear whether merging objects should be counted 
as one or two objects. Because we wish to have a sample 
which closely corresponds to the 'objects' seen in previous 
ground-based counts we smoothed the HST data with a l" 
FWHM Gaussian profile and then select objects above a 
fixed isophote of typically 25.5 mags/arcsec 2 , as if it was 
ground-based data. We then return to the unsmoothed data 
and classify each 'object' by detailed visual examination. 

Blue galaxies are seen in abundance at a mean redshift 
of ~ 0.5 where the WFPC2 pixel scale of 0.1" corresponds 
to ~ 0.4ft" 1 kpc. Since the PSF FWHM ~ 0.04" the data 
is under sampled and the effective PSF is close to a delta 
function in one pixel. Objects of galactic size are typically 
> 10 resolution elements across; the basic star-galaxy sep- 
aration is straight forward. Our total galaxy counts are 
shown in Figure 1(a) and compared with those from the 
ground (Tyson 1988, Lilly et al. 1991, Driver 1994). We use 
a small colour correction of /(F814W) = /(Cousins) — 0.07 
which represents the transformation of Harris et al. 1991 for 
a mean galaxy colour over 19 < / < 22 of (V — I) c = 1.5 
(Driver 1994, our data); we independently measured the 
likely I$i4 — I c colours with our galaxy spectral templates 



(described below) and got similar numbers. Figure 1(a) 
shows our counts are complete to / = 23 and agree well 
with the previous work. 

For our classified sample we consider all galaxies to 
I = 22 (301 total in all fields) irrespective of size; this allows 
us to obtain a clean magnitude limited sample which can be 
modeled in detail. At / = 22 (~= B = 24, K = 19) the 
excess over the no-evolution prediction is x2-4 (depending 
on local normalisation) and the median redshift for repre- 
sentative samples is ~ 0.5 (Lilly 1993). Obviously the more 
detailed a classification is attempted the more unreliable it 
becomes; differences between close types (e.g. Sa/Sb) be- 
come subjective. In this paper we are interested in only a 
basic classification of the fundamental morphological types, 
viz: 

i) Elliptical/SO galaxies. These show smooth profiles with 
no spiral features or disks. We see a range of sizes in- 
cluding a class of compact objects, contributing ~ 20% 
of the objects in each magnitude bin. 

ii) Spiral galaxies. We require these to show an extended 
disk with evident spiral arms, or to be clearly edge-on 
disks. Note in some cases the arms are seen to be dis- 
torted or show strong starbursts, but if the underlying 
structure is clearly a classic spiral we assign this class. 

iii) Irregular/peculiar galaxies. This category contains 
sources with peculiar or disturbed irregular structure. 
Many have low surface brightnesses compared to the 
previous two classes. We also include objects that ap- 
pear to be merging or closely-interacting (at sub-l" res- 
olution) . 

This first-order classification is fairly unambiguous and 
allows robust conclusions to be drawn. An extensive atlas 
of typical images of such MDS galaxies of various morpholo- 
gies has already been presented in Griffiths et al. 1994B. In 
this paper we now consider in detail the number-magnitude 
counts as a function of morphology, plotted in Figure l(b- 
d), to / < 22 and compare with no-evolution expectations. 

For our local luminosity function we use a determina- 
tion of Schechter parameters (M* , a) for early and late type 
galaxies taken from Loveday et al. (1992A) based on 1658 
APM galaxies to bj < 17. Zucca, Pozzetti & Zamorani 
(1994) have reanalysed this dataset under different assump- 
tions for classification errors at faint magnitudes and have 
suggested that the early type galaxies might have a flat faint 
end to their luminosity function rather than a turnover. We 
will consider this possibility. For a more detailed analy- 
sis we also use the morphological mix to bj < 16.7 given 
by Shanks et al. 1984, adjusting the weight of each type in 
the luminosity function to produce the observed flux-limited 
mix. Our K-corrections are based on spectral templates of 
Rocca-Volmerange & Guiderdoni (1988) which faithfully re- 
produce local galaxy UV-optical colours. 

The resulting models are plotted in Figure 1 for 
both normalisations: </>* = 0.015 h 3 Mpc -3 and </>* = 
0.03ft 3 Mpc -3 . For the ellipticals we also show the effect 
of a flat faint end to the luminosity function. Several points 
are of interest. Our total counts reproduce the excess noted 
in earlier work. However our counts for spiral and ellipti- 
cals show good agreement with the no-evolution prediction 
if we accept the higher normalisation. Moreover the slopes 
match the data indicating that any evolution would have to 
be mild. The effect of a flat faint end slope in the ellipti- 
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cal luminosity function affects only the very faint data and 
cannot yet be tested thoroughly. 

The most significant result is the steep rise in the num- 
ber of irregular/peculiar galaxies compared with their inci- 
dence in the local data. It is the increase which appears to 
produce the faint blue galaxy excess — ~ half the galaxy 
population at / = 22. 30% of them are either interacting or 
have multiple cores; these have the same steep slope as the 
isolated systems, confirming the tentative results of Colless 
et al. at a much higher resolution. 

How reliable are these results? We note that our clas- 
sification is independent of other groups in the MDS team, 
in particular Driver et al. (1994) have independently classi- 
fied a smaller sample of MDS galaxies and found a similar 
increasing proportion of irregular/peculiar galaxies at faint 
magnitudes. This suggests that the result is robust with 
respect to observer-observer classification errors. We also 
wished to assess the objective reliability; our most impor- 
tant concern, in view of our results, was that spiral disks 
might be hard to identify at high-redshifts leading to mis- 
classification as irregulars (due to spiral structure) or ellip- 
ticals (if not seen at all). Since for an object at z = 0.7 the 
2<r/pixel surface brightness limit in the HST data scales to 
typically 22.6 B mags/arcsec 2 at z = disks should still be 
easily visible. To check this we took a sample of CCD images 
of nearby galaxies and simulated observing the galaxies at 
z = 0.7 with WFPC2. A blind classification recovered the 
initial Hubble types; we generally found typical spiral disks 
were indeed very evident in the degraded data as expected. 
We conclude that there are no significant selection effects in 
our morphological classification due to the high redshift and 
faint magnitudes. 

The physical implication for the steep slope of the op- 
tical counts is clear. At bright magnitudes (B < 19) there 
must be some incompleteness in the Schmidt plate counts 
of elliptical and spiral galaxies to explain the apparent fast 
evolution claimed by Maddox et al. (1990). We note in- 
creasing support for high normalisation elsewhere: Glaze- 
brook et al. (1995B) find a local </>* = 0.026 ft 3 Mpc -3 from 
a near-infrared selected sample; Steidel et al. (1994) find 
</>* = 0.03 ft 3 Mpc -3 from a sample of luminous galaxies se- 
lected by their QSO absorption cross-sections. We feel that 
large-scale structures are unlikely to be the cause of this 
anomaly: it corresponds to physical scales of ~ 170 ft -1 Mpc 
and typical clustering variances a 2 are < 0.05 for scales 
> 60ft" 1 Mpc (e.g. Loveday et al. 1992B, Peacock 1991). 
Much more likely are explanations due to the high surface- 
brightness cut inherent in plate data (e.g. McGaugh 1994, 
Metcalfe et al. 1994B) or uncertainties in the local luminos- 
ity function (e.g. Koo et al.. 1993) or combinations thereof. 

However the rapidly evolving component is clearly iden- 
tified to B = 24 as being associated with objects of unusual 
morphology; and this population must contribute signifi- 
cantly to the large numbers of blue galaxies found at B > 24. 
Luminosity function analyses (Ellis et al. 1994) indicate 
that the evolving component at high redshift is luminous: 
Mb = —20 at z ~ 0.5, comparable to the elliptical/spiral 
populations. Evolution to lower luminosity, surface bright- 
ness or density (via mergers) must occur to remove these ob- 
jects by the present day. While a self-terminating starburst 
(Babul & Rees 1992) is one possibility; the high incidence of 
multiply-peaked systems suggests that mergers play a strong 



role. If so this must have little effect on the properties of 
the regular elliptical and spiral population (c.f. Dalcanton 
1993). The important discovery, from HST^s unique resolv- 
ing power, is the unambiguous identification of the evolving 
component to B = 24. Extensive spectroscopy is required 
to determine the detailed nature and evolutionary path of 
this population. 
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FIGURE CAPTIONS 



Figure 1. (a) The number-magnitude counts for all 
the MDS galaxies compared with the ground-based counts 
of Tyson (1988), Lilly et al. (1991) and Driver (1994). 
(b-d) The MDS galaxies split according to morphologi- 
cal classes. The lines are the no-evolution predictions de- 
scribed in the text for <f)* = 0.03 h 3 Mpc -3 (solid) and 
<f)* = 0.015 h 3 Mpc -3 (short dash) normalisations. The ex- 
tra long dash line on (b) shows the effect of a flat faint-end 
luminosity function on the E/S0 prediction for the high <f>* . 
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Fig. 1 
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